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i-Sh ■ Abstract. 



The Calogero-Moser systems are a series of interacting particle systems on one 
dimension that are both classically and quantum-mechanically integrable. Their 
integrability has been established through the use of Dunkl operators (a series of 
differential-difference operators that depend on the choice of an abstract set of vectors, 
or root system). At the same time, Dunkl operators are used to define a family of 
stochastic processes called Dunkl processes. We showed in a previous paper that when 
the coupling constant of interaction of the symmetric Dunkl process on the root system 
Aw-i goes to infinity (the freezing regime), its final configuration is proportional to the 
roots of the Hermite polynomials. It is also known that the positions of the particles 
of the Calogero-Moser system with particle exchange become fixed at the roots of the 
Hermite polynomials in the freezing regime. Although both systems present a freezing 
behaviour that depends on the roots of the Hermite polynomials, the reason for this 
similarity has been an open problem until now. In the present work, we introduce 
a new type of similarity transformation called the diffusion-scaling transformation, in 
which a new space variable is given by a diffusion-scaling variable constructed using the 
original space and time variables. We prove that the abstract Calogero-Moser system 
on an arbitrary root system is a diffusion-scaling transform of the Dunkl process on 
the same root system. With this, we prove that the similar freezing behaviour of the 
two systems on Ajy-i stems from their similar mathematical structure. 
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1. Introduction 

The Calogero-Moser (CM) systems are a family of exactly-solvable, interacting multiple 
particle one-dimensional systems. They have been of great interest in both physics and 
mathematics because they are one of the few interacting particle system families that 
can be solved exactly. The first such system was formulated by Calogero [JJ, and it 
is composed of N particles on a line confined by a harmonic potential, in which the 
particles repel each other through a potential proportional to the inverse of the squared 
distance between them. A similar system on the unit circle instead of the line was 
considered by Sutherland j2]. CM systems defined on the unit circle are commonly 
known as Calogero-Moser-Sutherland (CMS) systems. Later, Moser [3] established the 
integrability of the CM and CMS systems. These systems are categorized by the root 
system associated to them. A root system is a finite set of vectors that is closed under 
reflections; in other words, when one of its elements is reflected by another element, the 
result is an element of the root system as well (see Section [2]). A^-i is one such set 
of vectors. Olshanetsky and Perelomov proved that the systems defined in [1] and [2] 
are intimately related to the root system A^r-i, and that their analogues on all classical 
root systems are integrable [HE]. Polychronakos considered the inclusion of a particle 
exchange operator on the CM system first considered by Calogero [6]; this operator is 
also related to Ajv-i, so the CM system with particle exchange interaction is an extension 
of the CM system defined on pQ. He also studied its behavior when the potential 
strength parameter tends to infinity. This operation is known as the "freezing trick" 
[TJ. When the freezing trick is applied, the positions of the particles of the CM system 
with particle exchange interaction are fixed to the roots of the Hermite polynomials, so 
only the exchange interaction between particles remains as the non-constant part of the 
original system Hamiltonian. Because the particle exchange operator can be regarded as 
a spin interaction operator [8], the resulting system is called the Polychronakos- Frahm 
spin chain. Its integrability was consequently proved by Frahm [9j. The spectrum of 
the CM system with exchange operators on other root systems has been calculated by 
Khastgir, Pocklington and Sasaki [TO] . 

The integrability of the CM systems has been established by various methods, one of 
which makes use of Dunkl operators [UJ, [12]. These are differential difference operators 
that depend on the choice of a root system and a set of parameters called "multiplicites" , 
and they were defined for the study of multivariate orthogonal polynomials [TBI E]. 
When a similarity transformation is applied on the CM system Hamiltonian using its 
ground-state eigenfunction (see Section [3]), the result is an operator that has a very 
simple form when it is expressed in terms of the Dunkl operators [15]. At the same 
time, Dunkl operators are used to define Dunkl processes as a generalization of Brownian 
motion in the following sense: the Kolmogorov backward equation (KBE) that generates 
Brownian motion is the heat equation; if all the spatial derivatives in it are replaced 
by Dunkl operators, one obtains the Dunkl heat equation. Then, Dunkl processes are 
defined as the stochastic processes generated by the Dunkl heat equation [161 [17]. Dunkl 
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processes are discontinuous stochastic processes with a drift term in general, but their 
"radial part", or the part that is invariant under the reflection operators of the root 
system, is continuous. 

Demni [18] noted that the radial Dunkl process of type A N _i is equivalent to 
Dyson's Brownian motion [19]. We used this equivalence in our previous work [20] 
to give a physical interpretation to the Dunkl process of type It represents 

a system of N Brownian particles in one dimension that repel each other through a 
logarithmic potential and exchange positions spontaneously. Additionally, we found that 
the coupling constant for this interaction can be understood as the inverse temperature. 
More importantly, we used the mathematical tools of Dunkl processes to investigate 
the freezing properties of the radial Dunkl processes of type A N _i, and the equivalent 
Dyson's Brownian motion. We found that the (suitably scaled) trajectories of the 
Brownian particles are given by the roots of the iV-th Hermite polynomial multiplied by 
the square root of the process time when the temperature tends to zero. This is similar 
to the behaviour of the CM system with particle exchange interaction when the freezing 
trick is applied. The close relationship that exists between Dunkl operators and the CM 
systems indicates that the similarity in the freezing behaviour of the CM system and 
the radial Dunkl process of type A^-i is not accidental. The purpose of this work is 
to find a clear relationship between the CM systems and Dunkl processes that explains 
their similarity in the freezing regime. 

Our main result (Theorem [1]) is that one can use the diffusion- scaling transformation 
on the Kolmogorov forward equation (KFE) of any Dunkl process to obtain the 
Schrodinger equation of a CM system evolving on imaginary time, both on the same root 
system. The diffusion- scaling transformation is a similarity transformation in which a 
new space variable is defined by dividing the original space variable by the square root 
of the process time. Due to its importance in the derivation of our main result, we 
explain the diffusion-scaling transformation in more concrete terms as follows. Suppose 
that the frequency of the harmonic potential of the CM system is u. Then, we apply 
the variable substitution given by 



to the KFE of the Dunkl process. After that, we use the function e W ( T >£} } with 



to perform a similarity transformation on the result. Here, Wk(x) is a weight function 
that depends on the root system being considered, and its definition is given by flS]). 
Then, the result is that by applying the diffusion-scaling transformation to a Dunkl 
process (with space and time variables x and t, respectively) we obtain a CM system 
on the same root system (with variables £ and r). 

In Section [5J we review the definition of the reflection operator, the concept of root 
system and the definition of the multiplicity function. We also recall the Hamiltonian of 
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the CM system with exchange operators on an arbitrary root system R and the particular 
case of the CM system on Aat-i, as well as the definition of the Dunkl processes and their 
KBEs and KFEs. In Section El we review how a similarity transformation of the CM 
Hamiltonian using its ground-state eigenfunction yields a Dunkl Laplacian (a Laplacian 
operator with Dunkl operators instead of partial derivatives) plus a confinement term. 
In Section HI we review the main idea that gives the correspondence between Brownian 
motion and the quantum harmonic oscillator in one dimension. This idea, applied to the 
Dunkl process on a root system R, is the basis for the diffusion- scaling transformation 
we use to derive our main result, Theorem [1] in Section [51 That is, we prove that by 
using the diffusion-scaling transformation defined by the relations given in (TjQ) and (T5J), 
the KFE of a Dunkl process is transformed into the Schrodinger equation of the CM 
system. We explain the similar freezing behaviour of the Dunkl process and the CM 
system on the root system An_\ by applying Theorem CD in Section El We conclude by 
discussing these results and by outlining possible topics for future study in Section [71 

2. Preliminaries 

2.1. Root Systems and the Multiplicity Function 

We begin by defining the reflection operator by the equation 

x ■ ex 

a a x = x — 2 — — ex. (3) 

a 

This operator acts on the vector x by reflecting it through the hyperplane defined by the 
vector ex in N dimensions (ex, x e R ). Note that we write vectors with bold symbols, 
while we denote their squared norm by \ex\ 2 = ex ■ ex = a 2 . 

A root system is defined as a finite set of vectors such that, when its elements, 
called roots, are reflected by any root, the resulting vector also belongs to the set. In 
other words, a root system R is defined by the property that a a £ 6 R for any a, £ e R. 
In addition, a root system is called reduced if, for all ex e R, the statement hex e R 
with 6 6 R implies that b = ±1. We will assume that all the root systems considered 
here are reduced. We also assume that no root is the zero vector, because a reflection 
along the zero vector is undetermined. 

For every root system, there exists a set of base vectors with which every root is 
given by a linear combination of the base vectors with all the coefficients being either 
positive or negative. The roots in this base are called simple roots. This base is not 
unique, but it always divides the root system into two subsets: the positive subsystem 
R + (roots generated by positive coefficients) and the negative subsystem R_ (negative 
coefficients), and both subsystems have the same number of elements. 

An example of this is given by the root system B 2 = {±ei, ±e 2 , ±(e x — e 2 ), ±(e x + 
e 2 )} in M 2 , where denotes the i-th canonical base vector, % = 1, 2. If we choose e.\ and 
e 2 — ei as the simple roots, the positive subsystem becomes {ei, e 2 , e 2 + e\, e 2 — e{\. 
The negative subsystem is formed by the same roots multiplied by —1. In general, 
i2_ = {—ex : ex G (this is sometimes written as R = —R + .) 
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The multiplicity function, k(a), a G R, is a series of parameters, called 
multiplicities, assigned to each disjoint part of a root system in the following sense: 
if there exist roots ct,£, £ such that cr a £ = £, then k(£) = £;(£)• Consequently, the 
multiplicities assigned to two different roots is different only if they cannot be related 
by a series of reflections. 

In the case R = B 2 , there is no root that reflects the roots ±ei and ±e 2 into the 
roots ±(ei — e 2 ) and ±(ei+e2). Therefore, the multiplicity function for B 2 can take only 
two different values: k(±ei) = k(±e 2 ) = k% and k(±(e\ — e 2 )) = k(±(ei + 62)) = k 2 . 
That is, the root system B 2 has two multiplicities. In general, we will assume that all 
multiplicities are non-negative real numbers. 

It is useful to define the sum of multiplicities over the positive subsystem as 

7 = K<*) = \Y, k ( a )- ( 4 ) 

aeR + aeR 

Note that this value does not depend on the choice of simple subsystem. For example, 
for B 2 , 7 = 2{ki + k 2 ) for all possible positive subsystems. 

We will focus on the root system A^-\ in some of the later sections. This root 
system is given by A^-i = {(e.; — ej)/\/2 : 1 < i 7^ j < N}. Choosing the simple roots 

(Ci+x - ei)/V2, i = l,...,JV-l (5) 

produces the positive subsystem 

Ajv_ li+ = {onj = (e, t - ej)/V2 :l<j<i<N}. (6) 

Additionally, this root system has only one multiplicity, k, and we have 
iV 

7= -(AT -i) k . (7) 

Most importantly, the effect of the reflections along the roots of A^-i is that of 
exchanging the reflected vector's components. Using the notation = o~ aiP the 
expression a^x denotes the vector obtained by exchanging the z-th and j-th components 
of the vector x. For details on the above, see Appendix A on [2U] . 
Finally, we define the following weight function: 

w k (x) = Yl \a-x\ k ( a \ (8) 

This function is invariant under reflections along the roots of R, and it is harmonic, i.e., 

N f) 2 

A^w k (x) = Y,g^Mx) = 0. (9) 

i=l % 

This function, w k (x), will be used in the transformation considered in Theorem [1] in 
Section 
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2.2. Calogero-Moser systems 

Under an arbitrary root system R, the CM systems on a line with a harmonic background 
potential and an inverse-square repulsion potential are given by the Hamiltonian (see, 
e.g., HO]) 

njR l/vOr) i y a 2 k{a)[k(a) - a a ] a; 2 A 2 

aGi?+ V ' i=l 

where all particles have been chosen to be of unit mass, and we have taken H = 1. In 
particular, when the root system is Atv_i and the harmonic frequency u is set equal to 
the multiplicity k, the CM system with particle exchange interaction considered in [7] 
is recovered: 

= £ + (ii) 



l<i<j<N y 1 JJ i=l 



2. 3. Dunkl processes 



We define the gradient operator as = (^-, • • • , af - )- Under an arbitrary root 

system R, the Dunkl operator in the direction £ is given by [13J 

TJ(x) = t ■ VW/(«) + J2 k{<*) f{X) ~ f{GaX) ct ■ & (12) 

a&R+ 

that is, it consists of a combination of a directional derivative and a sum of difference 
terms taken along the reflections given by the roots in the positive subsystem. If we 
choose £ = e», we use the notation T Ei = Ti, and we write 

T if (x) = *-f{x) + £ W S{x) i n ° aX) a<. (13) 
The Dunkl Laplacian is given by [2] 

(14) 



z — ' L — ' a x z — ' (a • xY 

1=1 o:ei?+ aei?.+ v 7 

and with it, the Dunkl heat equation is defined by [16] 



N 



^u(t,x) = ~J2Tlu(t,x), (15) 

i=l 

where u(t, x) is a sufficiently well-behaved function. Dunkl processes are defined as the 
stochastic processes whose KBE is the Dunkl heat equation. That is, their transition 
probability density (TPD) pk(t,x\x') obeys the equation [141 [T7] 

-p k {t, x\x') = -A** )p fc (t, x\x') + 22 k ( a) 



a£R+ 

E a 2 p k (t,x\x') - p k (t,x\a a x') 
k ( a H / ~A2 > ( 16 ) 
2 (a • arr 



Calogero-Moser Systems as a Diffusion- Scaling Transform of Dunkl Processes 



7 



where we denote the initial condition using primed variables. Using the explicit form 
of p k (t, x\x') (see, e.g., [21]) and the KBE above, it can be shown directly that the 
corresponding KFE is given by 

d_ 
di 1 



;Pk (t,x\x') = -A^p k (fx\x') - Ha) ™ K J 1 ; 



ST u ^ 2 Pk{t,x\x')+p k {t,a a x\x') 

+ 2. k ^Y (^xj • (17) 

From the two Kolmogorov equations, one can read off the terms on the rhs, from left to 
right, as a diffusion term, a drift term and a spontaneous jump term, respectively. The 
drift and jump terms are clearly dependent on the root system. In particular, the Dunkl 
process of type A^-\ can be understood as the Brownian motion in one dimension of N 
particles which repel each other through a logarithmic potential and exchange positions 
spontaneously [20J. This interpretation is reached by writing down the corresponding 
KBE [H], 

9 i N N k o 

—p k {t, x\x') = -A^ Pk (t, x\x') + rf^d.~faJp k & x ^ 

i=l j=l: % 3 1 

k \^\^ gfcfe x \ x ') ~ gfcfo x \ a v x ') (-to) 

1=1 3=1: 3 »' 

and noting that this equation is similar to the KBE of Dyson's Brownian motion |22j . 

Finally, we describe radial Dunkl processes. These are the stochastic processes 
that result from symmetrizing the initial condition of p k (t, x\x') under the reflections 
generated by the root system. In this case, the jump term on the KBE vanishes, and 
the corresponding Dunkl process becomes a continuous stochastic processes. As noted 
in [IS], the radial Dunkl process of type A^-i corresponds to Dyson's Brownian motion, 
with (3 = 2k. Here, (3 is Dyson's beta parameter, and it is equal to 1,2 and 4 for the 
Gaussian orthogonal, unitary and symplectic random matrix ensembles, respectively 

[121 123]. 



3. Similarity Transformation of the Calogero-Moser Systems 

Dunkl operators have been used as a tool to prove the integrability of the CM systems 
[T2] . More specifically, it has been shown under several root systems [16] that after 
applying a similarity transformation (using the ground state eigenfunction), the CM 
system Hamiltonian is expressed as a Dunkl Laplacian minus a term of the form cc-V^. 
For R = An-i, given the function W(x) = J2iLi x l/2 ~ Z^i<j ^ n l x * — x j\-> ^ e operator 
equation 

- e™(n A c ^ - E^)e- kw = \Y, T "- k Y, x ^. (19) 

i=l 7=1 J 
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is obtained. Here, the ground-state energy is given by E^' 1 = [kN + k 2 N(N — l)]/2. 
One can then proceed to find polynomial eigenfunctions for the rhs of this expression 
as stated in [15] and shown in [24J. 

The equation above suggests that there must be a way to transform the CM 
Hamiltonian of type A^-i into an expression proportional to the Dunkl Laplacian of 
type Ajv-i- However, we would like to obtain a relation like f lT9|) without the second 
term on its rhs. 



4. Brownian Motion and the Quantum Harmonic Oscillator 

A similar relationship appears between one-dimensional Brownian motion and the 
quantum harmonic oscillator; consider the Hamiltonian for the harmonic oscillator with 
m = lo = H = 1, "Hh = ~\j^2 + \ x ' 1 ■ A quick calculation yields 

-e^ ( « H _ 1/2)e -^ = l^_^. (20) 

In other words, the harmonic oscillator can be transformed into a Laplacian minus the 
term x-^-. The stochastic process generated by the rhs of (1201) has the KBE 

d Id 2 d 

and it is called an Ornstein-Uhlenbeck process [25]; it represents a Brownian motion 
confined by a quadratic potential. Therefore, we can understand the rhs of f)19p as the 
generator of a Dunkl process in a harmonic potential. 

There also exists a way to transform the KBE of a Brownian motion (the heat 
equation, which is also its KFE) into a form that resembles the Schrodinger equation of 
the harmonic oscillator [22]. Consider a Brownian motion whose probability density is 
given by the function u(t,x) with initial condition u(0,x) = 5{x — x'), t > 0. Then, its 
KFE is 

d , . Id 2 
m U ^ x) = 2Tx 2 
Applying the variable substitution 



-u(t,x) = -—u(t,x). (22) 



r,C)=(lnt,4=) (23) 



<2t 

and the similarity transformation u{t, () = exp[— (r + ( 2 )/2]U(r, (), (|22|) becomes 

d -U(T,C) = -l(Hn-l)u(T,C), (24) 



Or y 2 V 2, 

with x replaced by ( in % H - This last equation shows that, under the transformation 
considered, Brownian motion corresponds to the harmonic oscillator evolving in 
imaginary time. 

Intuitively speaking, this transformation succeeds in relating Brownian motion and 
the harmonic oscillator because of the variable transformation (|23|) . Because Brownian 
motion is a free diffusion process, it is unbounded in space, whereas the harmonic 
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oscillator is not. Dimensional analysis on the heat equation indicates that x scales as 
y/i, so we expect this diffusion process to be bounded when such a scaling is applied on 
x. This scaling yields a process similar to the Ornstein-Uhlenbeck process described by 
(pTj) . After that, applying a suitable similarity transformation produces the harmonic 
oscillator Hamiltonian. 



5. Dunkl Processes and CM Systems on an Arbitrary Root System 

The ideas of the previous section are applicable to Dunkl processes as follows: consider 
now a straightforward generalization of the substitution (j23|) given by ([1]), that is, 
consider a new space variable defined by diffusion-scaling the original space variable. 
Because the scaling is isotropic, it should have the same effect of fl23l) when applied 
on the Dunkl process on R. Hence, the diffusion- scaling transformation should relate 
Dunkl processes and the CM systems in the same way that Brownian motion and the 
quantum harmonic oscillator are related, as seen in the previous section. That is, we 
expect to be able to transform the Dunkl process on R into the CM system on the 
same root system using the diffusion-scaling transformation, that is, by applying the 
similarity transformation 

u(rX) = ex P [-W(rX)]U(rX) (25) 

with W(t, C) given by ([2]) after the substitution (JT]). However, we must be careful with 
the type of Kolmogorov equation which we choose to transform. Brownian motion has 
the same equation for its KBE and KFE, but fl22|) represents its forward evolution, so 
we should transform the KFE ( Tl7l) . We prove that the diffusion-scaling transformation 
indeed transforms Dunkl processes into CM systems in the following theorem. 

Theorem 1. The diffusion-scaling transformation given by (QJ), (E)) and l[25\) transforms 
the Dunkl process on the root system R into the Calogero-Moser system with harmonic 
confinement on the same root system evolving in imaginary time. 

Proof. We begin by considering a Dunkl process whose distribution u(t, x) obeys both 
the initial condition 

N 

u{0,x) = 5 N {x- x') = Y[5{x i -x' i ), t>0, (26) 

i=l 

and the KFE (TT7|) . We write down the derivatives in time and space in terms of the 
new variables as follows: 

dt 2utdr 2t s 

— = — - (27) 

dXi y^JtdQ' 

Note that these derivative transformations are undefined for t = 0. Thus, we have chosen 
t > to guarantee the consistency of the transformation. We insert these derivatives in 
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(fT7|) to obtain 



t E <'''i t .t'' ,{l ^'^,fl. P8) 



The differential operators above are transformed by (125]) as follows: 

e = a uN > 



e w A«'e- w = A«) + 2 ( £ " • V (c) + a; 2 C 2 - (2 7 + N) 



ae-R.4 



We insert the above into ( 1281) and obtain 
d 



^ U(t, = Ia^U(t, C) + | [2 7 + N - u,C 2 ]tf(r, C) 

Ea 2 Ma) rr , 



- 2^ 2.— (a.c)(£-C) ( ' C) - ( } 

The double sum in the bottom term of the equation above can be simplified because all 
the terms where a/^ cancel among themselves. We prove this in the Appendix. We 
finally obtain 

^U(r, C) = \a^U(t, C) + I [2 7 + N - co( 2 ]U(t, 

or, denoting the ground-state energy by = lu('j-\-N/2) and using with £ instead 
of x, 

~ §^U(t, = [H% M - £ >(r, C), (32) 

as desired. □ 

Remark: this proof involves only straightforward calculations, with the notable 
exception of the step required to obtain (I3T]) . This is perhaps the most important part 
of the proof, and it is not trivial. In the following section we focus on the root system 
An_i and we provide an explicit example of the calculations above. In that case, the 
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simplification needed to obtain ( 131 j) is straightforward. However, the general case is 
more involved, and it is detailed in the Appendix. 

This theorem implies that CM systems and Dunkl processes share a common 
mathematical structure and therefore behave similarly. In particular, their freezing 
behaviour should be similar regardless of the root system under consideration. However, 
one must be careful when taking the freezing limit because there is more than one 
parameter involved, leading to multiple ways to take the limit (the multiplicity function 
k{ot) represents more than one parameter in general, and there is also the harmonic 
frequency u>.) This result also implies that the wealth of information available on the 
spectrum and eigenfunctions of the CM systems can be used to study Dunkl processes. 

Note that Theorem [1] only requires that oj > 0. If u = 0, there is no need to use 
the diffusion scaling (QQ), and one may simply apply a similarity transformation on the 
Dunkl process to obtain the unconfmed CM system on the same root system. 



6. Dunkl Process and Calogero-Moser System of Type An_i 

Having proved Theorem [H we may consider the particular case of the root system 
An^i as a corollary. In this case, we map the Dunkl process of type Aw-i to the CM 
system with particle exchange interaction described by "H^ivf 1 m dHJ)? so we apply the 
diffusion-scaling transformation using (0Q) with u = k and (|25l) with 

W(rX) = ^( 2 -k\n H iQ-Q+kNr (33) 

l<i<j<N 

plus a constant term that vanishes in the transformation. 

Corollary 2 (Particular case of Theorem [1]). The diffusion-scaling transformation given 
by (GP; ( E5P and (E3|) with uj = k transforms the Dunkl process of type An_\ into the 
Calogero-Moser system with particle exchange under harmonic confinement evolving in 
imaginary time. 

Proof. It suffices to follow the same procedure of Theorem [1] applied to An_\ to obtain 
the result. As before, we consider the distribution u(t,x) with initial condition ( 126|) . 
The KFE (JTZJ) becomes 

8 1 N N k d 

-) = 2 AW " (t - ^'-EE -) 

i=\ j = l; 1 J 1 

i_ t V 1 V 1 u(t,x) + u(t,aijx) 

+ 2^-^f ( Xj - Xif ' 1 ' 

i=l j=l: ^ J % > 

The expressions (J2ZD, ®, (J29D and §^ become 

dt 2ktdr 2t s 
did 
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d_ 



1 N N 



k d 



i=l j = l: 



Ci ~ Cj dG 



u(t,C) 



N N 



1=1 j = l: 



or or 



+ \ E E " (T '(^f j0 + < ■ V«'«(r, o, 



N 



^ Ea_£. 
j = l: ^ ^ 

AT r AT 



e »- A (0 e -H' 



i=l: ^ ^ 



50 



AT N 



kC-kN(N-l)-N + J2Yl 



k-l 



1=1 J 



and 



|-£/(t, C) = l^U(r, C) + -1)+N- kC 2 ]U(r, Q 



i=l j=l: 
N N N 



(o - ay 



t/(r, C) 
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(36) 



(37) 



(38) 



1= t^.fe (C< — OXCi — CO 

respectively. The triple sum in the bottom term is simplified by considering the case 
i / j / I, where we have the following three term sum: 

1 1 1 



(Q-Q(Q-Q) (0-0)(Ci-G) (O-CO(O-O) 
_ G - - Cj + + Ci - Cj _ n ( o q \ 
(C-0)(C-C)(Cj-0) ' 1 j 

Therefore, the only remaining terms are those in which i ^ j = I; this is a particular 
case of the simplification of the same term in (1301) . as proved in the Appendix. Then, 
(ED and 032]) become 

N N 



Or 



d U(r, O = \ ^U(r, C) - | E E ^'g:^ 

A; 2 u „ jfeJV, 



1=1 j = l: 



yC 2 ^, C) + '^HN - 1) + l]U(r, C) 



(40) 
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and 

- §pU(r, C) = - E^]U(r, C), (41) 

respectively, proving the statement. □ 

This corollary explains the similarity of the freezing behaviour of the radial Dunkl 
process of type Ajy-i and the CM system on Ajv-i- On one hand, the transition 
probability density of the radial Dunkl processes freezes to [20] 

lim p s k (t, \fkv\x')h N/2 = ?4(t, v)= V 6 N [v - V2tpz N ], (42) 

pes N 

where is the vector whose components are the (ordered) roots of the iV-th Hermite 
polynomial, defined by [26] 

H N (x) = (-l)^ 2 -^-* 2 ). (43) 

On the other hand, under the same limit, the CM system with particle exchange 
interaction freezes to the Polychronakos-Frahm (PF) spin chain [7]. That is, the 
positions of the particles freeze at the roots of the iV-th Hermite polynomial and the 
only remaining dynamical part of Hq 1 ^" 1 is the PF spin chain Hamiltonian, 

Kpf= £ (44) 

If we recall the diffusion-scaling substitution (JTJ), we realize that 

xv , . 

C = -j= = 45 

because x = ykv in (142 p . This equation is consistent with the freezing regimes of 
both systems (£ — > and v — > \[2tz^ when the freezing limit is taken), so the factor 
of y/2t is accounted for by the diffusion-scaling transformation. In other words, the 
similar behaviour of both systems in the freezing limit goes beyond simply multiplying 
or dividing a factor of v2i to the particles' frozen positions. It is a consequence of 
the fact that the CM system with particle exchange interaction is the diffusion-scaling 
transform of the Dunkl process of type Aw-i- Therefore, we can say that because 
one system is the diffusion-scaling transform of the other, they behave similarly at any 
temperature and in particular, they freeze in the same way. 



7. Concluding Remarks 

We have established a correspondence between Dunkl processes and CM systems with 
harmonic confinement on a line with the use of the diffusion- scaling transformation. This 
strategy works because the variable substitution confines the Dunkl process equally 
in all directions by producing a restoring drift term (the second term on the rhs of 
the first line of (j27|) ). Once this drift term is present, it only remains to perform a 
similar operation to the one given in f fl9|) to obtain the desired result. In our case, 
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that operation is a similarity transformation based on a function proportional to the 
ground-state eigenfunction of the CM system, exp[— W(r, £)] [TO] . 

By establishing this correspondence, we conclude that these two types of systems 
share similar characteristics. In particular, their freezing limits must be similar in the 
sense that, if the particles in a CM system freeze at the positions z, then the scaled final 
position vector of the corresponding radial Dunkl process is given by v = \/2iz, where 
we assume that u = k(a) = k and v = x/\/k before taking the limit k — > oo. The 
particular case of the CM system of type B N is studied in [27] where it is shown that 
when u = ki = k 2 = k, the system freezes at the roots of the Laguerre polynomials. 
Therefore, we expect that the radial Dunkl process of type Bn should freeze to a scaled 
final position proportional to the Laguerre roots. This is a calculation we plan to tackle 
in the near future. 

Concerning the circular case, T. Kimura has pointed out that the the freezing 
trick has also been applied to the CMS (or circular CM) systems. In that case, the 
freezing trick produces spin chains of Haldane-Shastry type [28, 29J. This prompts the 
possibility of constructing a mapping from Dunkl processes defined on the unit circle 
to CMS systems similar to the diffusion-scaling transformation discussed here. At first 
glance, we expect the mapping to be fairly straightforward, because both systems are 
defined on a space of finite size. Hence, no diffusion scaling would be needed and only 
a similarity transformation should be required. However, circular Dunkl processes have 
not yet been defined to the best of our knowledge. We suspect that the trigonometric 
Dunkl operators considered by Cherednik [30], Heckman [31] and Opdam [32] should 
be at the core of such a definition. This is a problem we leave open for future work. 

As noted before, there exist different ways of taking the freezing limit. So far, 
we have considered the cases where all parameters are locked to the same value and 
then brought to infinity. However, little is known about freezing limits where not all 
the multiplicities tend to infinity, and in consequence, their physical implications are 
unknown as well. They have been studied in the context of the generalized Bessel 
function on certain root systems [33, 34J. The generalized Bessel function is relevant 
in the sense that it is a part of the TPD of radial Dunkl processes [18} 135} 136] . It has 
been shown that this function undergoes a transition from the root system Bn with 
multiplicities ki, k 2 to the root system An_i with multiplicity k = k 2 when ki tends to 
infinity and its arguments are properly scaled. It is of interest to test whether this kind 
of transition occurs in the context of Dunkl processes and CM systems, and we plan to 
investigate this matter starting with the root system B^. 

Finally, and following the ideas in [22J, the fact that the CM systems with harmonic 
confinement are the diffusion- scaling transform of the Dunkl processes on the same 
root system allows us to use the quantum mechanics of the CM systems to study 
Dunkl processes. In particular, the harmonic confinement produces a discrete basis 
of eigenf unctions which can be used to represent the TPD of Dunkl processes. Hence, 
we believe that the existing results on the CM systems may provide further insight into 
the theory of Dunkl processes and their applications in physics. 
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Appendix 

We prove two lemmas necessary to complete the proof of Theorem [U The first lemma 
corresponds to Theorem 4.2.4 of [14] , and it involves a polynomial called the discriminant 
of a (reduced) root system R, defined by 

a R (x) = \{v x (A.l) 

v€R + 

Lemma 3. The discriminant of R is the minimum-order polynomial that obeys the 
alternating property 

a R (a u x) = -a R {x) (A. 2) 

for all u G R, up to a constant factor. 

Proof. By definition, a root system is invariant under reflection along any of its elements. 
Hence, we can divide R + into the following three sets: 

E% = {v eR + : a u v = v}, 

E% = {v G R + : V ^ v s.t. a u v = ±v'} (A.3) 

and {u}. That is, when a root of R + is reflected along u, it is either unchanged 
(orthogonal), it is reflected onto another root of R, which is in R + or in R_ = —R+, 
or it is turned into its negative, i.e., it is u itself. By definition, the roots in E^ obey 
the property that if v G E% , there is a root v' such that o~ u v = ±i/, so applying the 
reflection a u on this last equation yields 

o u v' = ±v. (A.4) 

That is, v' G , so there is an even number of roots in , and they form pairs such 
that 

[p u v ■ x)[o u v' ■ x) = [v' ■ x){v ■ x). (A.5) 
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Now, we can write 

o>r{(TuX) = } [ °~ u v ■ x = a u u ■ x ] [ o u v ■ x a u v ■ x 

— — u ■ x J | v ■ x J J v ■ x = — J [ v ■ x = —aft(x) (A. 6) 

The reason for the first equality is that a u x ■ y = x ■ o u y for any vectors x and y, with 
u / 0. This equation holds for all u G R. Moreover, if we were to take any one of 
the factors that make up a r (x) away, say the factor u' ■ x, it would no longer be an 
alternating polynomial along u' . This completes the proof. □ 

The following lemma is a particular case of Lemma 4.4.6 of [14J. It requires Lemma El 
and it completes Theorem [IJ 



Lemma 4. The double sum in ptyl is given by 

y- £ a . s JMMfL = £ „ 2 Jie^. (A7) 

^ ^ (a-x)(£-x) ^ (a ■ x) 2 

Proof. Let us consider two different roots of R + , say, u and v. We denote by w uv the 
rotation obtained by composing the reflections along u and v, that is, w uv = cr u a v . Let 
us also define the function 

/<-->- Jji^f- 

We will prove that f(x,w uv ) = 0. For this purpose, we define the function 

{x)f(x ). (A.9) 

Here, Span(u, v) represents the vector space (plane) generated by the pair of vectors 
{u,v}. Let us calculate f(a z x,w uv ), with z G R D Span(w,t;), that is, z is a root 
contained in the plane of the rotation w uv . 

^ (a ■ £)k(a)k(£) ^ (a z a ■ a z ^)k(a z a)k(a z ^) 

J{a z x,w uv ) = > — -= > — 

' (a ■ o z x £ • o z x) (a z a. ■ x)(a z t ■ x) 

On the second line, we have used the variable substitution a z ct = a', and the analogue 
for A direct calculation reveals that a ai!Ct i = o- z a a /a z , so 

o~a z oL'0-a^> = o- z a a >a z a z a^a z = cr z a a ra^a z , (A. 11) 

and the condition o~ az a'0~a z £' = w uv becomes 

o~(io~ a i = w uv . (A. 12) 
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Hence, 



f(a z x,w. 



uv 



)= £ 




ot' ■ i'k{a!)k{i') 
(a' • x)(£' ■ x) 



f(x,w uv ) 



(A.13) 



because each term in the sum is unchanged when and ex.' are exchanged. Now, we 
know that 



because of Lemma [31 so g(x,w uv ) also has the alternating property. However, 
a RnSp&n(u,v)(x) is the alternating polynomial on R fl Span(tt, v) of minimum degree 
(|.R + nSpan(it, v)\), and g(x, w uv ) is a polynomial of smaller degree, |i? + flSpan(it, v)\ — 
2. This contradicts Lemma |3] unless g(x,w uv ) is identically equal to zero. Therefore, 
f(x,w uv ) must also be equal to zero. 

Repeating the argument above for all possible pairs of different u,v e R + , one can 
show that the terms where a 7^ £ in ( 1A.7I) vanish, and only the terms with a. = £ remain. 
The remaining terms do not vanish because, in the argument above, if v = u then 
a mspa.n(u,v){. x ) is a linear function of x. Hence, g(x,w uv ) is no longer a polynomial, so 
we cannot use Lemma[3]to arrive to a contradiction. With this, the proof is complete. □ 
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